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Abstract 

A new method of interference contributions accounting is proposed. This method makes it possible 
to solve the problem of accounting the particles identity in the final state of the scattering process. This 
problem requires the calculation of a huge number of multidimensional integrals that express the scattering 
cross-section through the square modulus of the scattering amplitude. Our method is based on the physical 
features of the scattering amplitude. We show that the set of all interference contributions may be divided 
into subsets with the similar contributions. Therefore one does not need to calculate each of them separately. 
Within this approach it is possible to describe inelastic scattering processes with large number of secondary 
particles. The method was also verified within the model using the Laplace method for the calculation of 
multidimensional integrals in order to compute the inclusive rapidity cross-sections. A qualitative description 
of experimental data was obtained. 

An important result of this work is the explanation of peaks behavior in the inclusive rapidity distributions 
and their possible interference origin. 


1 Introduction 

There are two problems when considering an inelastic scattering processes: the first problem is the calculation 
of a multidimensional integral (1); and the second problem is the large number of such integrals. 
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The first problem may be solved using the Laplace method. 

Now let us consider the so called disconnected diagram (Fig. 1). In the language of diagrams: we need to 
connect the lines of secondary particles pi,p 2 , ■■■,n to the diagram vertices l,2,...,n in all possible ways. For 
each way of connection we need to calculate the interference contribution i.e. the integral (1). There are n\ 
possible ways of such connections. It means that we need to calculate n\ multidimentional integrals. It is not a 
problem for a small number of secondary particles (e.g. for five or seven). However, for the process with a large 
number (20,50,60...) of secondary particles, it is impossible to calculate all interference contributions directly. 


Secondary particles number 

Number of interference contributions 

n = 5 

n! = 120 

n = 7 

n! = 5040 

n = 15 

n! = 13 • IQii 

n = 60 

n! = 8.32 • 10®i 


Table 1: The number of the interference contributions for the different secondary particles number. 

Obviously, it is impossible to build all the permutations even for 30 particles nowadays, and thus to ac¬ 
count all 30! possible interference contributions. Our method makes it possible to calculate the interference 
contributions for up to hundred secondary particles. 


1 











0 


1 

2 

3 

4 

5 


Pi 

P2 

P3 

P4 

Ps 



Figure 1: Disconnected Feynman diagram. 


2 Laplace method 

The scattering amplitude A (Pi, P 2 , P 3 , P 4 ,Pi,P 2 , ■ ■ ■ ,Pn) depends on the four-momenta of the secondary parit- 
cles. It is convenient to introduce the (3n-|- 2)-dimentional vector X, which includes all the necessary variables 

yi : 1/2 5 ■ ■ - ; Vn 1 Plxt P2x ; ■ ■ • ; Pnx 7 Ply 7 P2y 7 7 Pny and Pax — (P^x Pax^ 7 Pay — (P^y P^y^ l‘^- And SO let US 

write A(y/s, X). 

The scattering amplitude may be given in the following form; 

A(v^,X) = (2) 

As it was shown in [1], the scattering amplitude A(^,X) has a maximum point Xq at fixed energy ^/s, 
and by applying Taylor expansion near the maximum point Xq up to second derivatives: 

A(v^, X) = A(yG, (3) 

where 


Dab = 


d‘^lnA(i/s, X) 


dXadXt 


X=Xo 


(4) 


The diagonalization of matrix Dab yields several one-dimentional integrals instead of a single multidimen- 
tional one. 


3 Definitions and important details 

First of all we need to introduce the definition of the ’’normal connection” or ’’normal order”. Let us connect 
first line to the first vertex, second line to the second vertex and so on. From now on let us call it a ’’normal 
order” or a ’’normal connection”. Now we can consider all the possible connections as the permutations of the 
seconrady particles in normal order. 

It is important that the secondary particle rapidities make up the arithmetic progression at the maximum 
point of scattering amplitude (Fig. 3 right), while Umax = —Umin and the central particle has a zero rapidity. 
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Figure 2: Normal connection of secondary particles (left), and secondary particles rapidities at the maxtimum 
point (right). 


4 The essence of the new method 


The permutations of the secondary particles (different connections) at low energies lead to the same maximum 
point of scattering amplitude (Fig. 3), thus there are many similar contributions which are easy to calculate. 
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Figure 3: Low energies: different connections but the same maximum point. 

This is not true for the high energies. However, at high energies it is possible to combine the secondary 
particles into the groups in such a way that the permutations of the particles inside the groups do not change 
the maximum point. It means, that the rapidities of the particles from the same group are equal. 

The problem of accounting of all n! interference contributions may be solved in two steps: 

1. The first step is the consideration of all possible permutations inside the groups - inner permutations 

2. The second step is the consideration of all possible ways (methods) of groups filling - external permutations. 


A{y/s,P,p) = E E a{'/s,P3,P4,Pl,P2,--,Pn) (5) 

external perm, inner perm. 

From now on the Laplace method may be applied to the partial sum of the interference contributions: 

^ a{y/s,P3,P4,Pl,P2,-,Pn) (6) 

inner perm. 

All inner permutations may be accounted by the combinatorial coefficients. Each external permutation may 
be represented by x A-dimentional matrix. 
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M,k (7) 

The matrix element Mik shows how many patricles are taken from the group of the normal order and 
then pushed to the group of the new permutation. Each of these contributions has a weighting factor Pm '■ 


ft'=nnc, 

i=li=l 


Mi. 




' Mil 


( 8 ) 


In this way the number of necessary calculations of the multidimentional integrals is reduced by finding 
the similar interference contributions and accounting all of them by combinatorial coefficients. So it becomes 
possible to calculate all the interference contributions for the processes with a large number (up to hundred) of 
the secondary particles. 


5 Results and Conclusions 


As a result of this method application, a qualitative description of inclusive rapidity distribution was achieved. 



Figure 4: Inclusive rapidity distribution. The experimental data are taken from [2]. 

It is important to note that there is a single peak for low energies (Fig. 4), which splits into two peaks more 
and more as energy grows. 

5.1 An explanation of peaks behavionr 

Let us consider the so called cut diagram (Fig. 5). There are two types of lines connecting the secondary particle 
lines from the left side to the right side: 

• Green (dotted) lines connect the particle lines from the central group on left side to the particles from the 
central group on the right side (Fig. 5). It is these lines that produce the maximum at zero point, since 
the particles from the central groups have zero rapidities. 

• Red lines connect the particle lines from/to the non-central groups (Fig. 5). The lines of this type shift 
the maximum from zero, because the particles in non-central groups have nonzero rapidities. 

There are two factors impact on the peaks behaviour with the energy y/s growth: 

• The possible number of the secondary particles (n) increases, and thus the number of the red and green 
lines increases. 

• A difference of the rapidities arithmetic progression increases. 

However, the number of red lines grows faster than the number of green lines as the secondary particles 
number increases. This fact provides the peak splitting. 
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